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ABSTRACT 

The discovery of the second highly r-process-enhanced, extremely metal-poor star, CS 31082- 
001 ([Fe/H] = —2.9) has provided a powerful new tool for age determination, by virtue of the 
detection and measurement of the radioactive species uranium and thorium. Because the half-life 
of 238 U is one-third that of 232 Th, the U-Th pair can, in principle, provide a far more precise 
cosmochronometer than the Th-Eu pair that has been used in previous investigations. In the 
application of this chronometer, the age of (the progenitor of) CS 31082-001 can be regarded as 
the minimum age of the Galaxy, and hence of the universe. 

One of the serious limitations of this approach, however, is that predictions of the production 
ratio of U and Th have not been made in the context of a realistic astrophysical model of the 
r-process. We have endeavored to produce such a model, based on the "neutrino winds" that are 
expected to arise from the nascent neutron star of a core-collapse supernova. In this model, the 
proto-neutron star mass and the (asymptotic) neutrino sphere radius are assumed to be 2.0Af Q 
and 10 km, respectively. Recent hydrodynamic studies indicate that there may exist difficulties 
in obtaining such a compact (massive and/or small in radius) remnant. Nevertheless, we utilize 
this set of parameter choices since previous work suggests that the third r-process peak (and 
thus U and Th) is hardly reached when one adopts a less compact proto-neutron star in the 
framework of the neutrino-wind scenario. The temperature and density histories of the material 
involved in the neutron-capture processes are obtained with the assumption of a steady flow of 
the neutrino-powered winds, with general relativistic effects taken into account. The electron 
fraction is taken to be a free parameter, constant with time. The r-process nucleosynthesis in 
these trajectories is calculated with a nuclear reaction network code including actinides up to 
Z = 100. 

The mass-integrated r-process yields, obtained by assuming a simple time evolution of the 
neutrino luminosity, are compared to the available spectroscopic elemental abundance data of 
CS 31082-001. As a result, the "age" of this star is determined to be 14.1 ± 2.5 Gyr, in excellent 
agreement with lower limits on the age of the universe estimated by other dating techniques, 
as well as with other stellar radioactive age estimates. Future measurements of Pt and Pb in 
this star, as well as expansion of searches for additional r-process-enhanced, metal-poor stars 
(especially those in which both U and Th are measurable), are of special importance to constrain 
the current astrophysical models for the r-process. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — supernovae: 
general 
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1. INTRODUCTION 

Uranium and thorium are regarded as poten- 
tially useful cosmochronometers because their long 
radioactive decay half- lives ( 232 Th: 14.05 Gyr, 
238 U: 4.468 Gyr) arc significant fractions of the 
expected age of the universe. All of the actinides, 
including uranium and thorium, are thought to 
be pure r-process elements, since the alternative 
neutron-capture process, the s-process, terminates 
at 209 Bi. The excellent agreement of the relative 
abundances of neutron-capture elements in the ex- 
tremely metal-poor ([Fe/H] = —3.1), r-process- 
element-enhanced star CS 22892-052 with the so- 
lar r-process pattern initially suggested that tho- 
rium might serve as a precise cosmochronometer 
(Sneden et al. 1996). The time that has passed 
since the production of the thorium that is now 
observed in the outer atmosphere of such an old 
halo star can be regarded as a hard lower limit on 
the age of the universe. 1 Sneden et al. (1996) es- 
timated the age of this star to be 15.2 ± 3.7 Gyr, 
comparing the measured Th/Eu ratio with an ini- 
tial production ratio constrained by theory. This 
age limit, albeit with a rather large error bar, is 
in good agreement with the age limits derived by 
completely different techniques, e.g., for globular 
clusters (12.9 ± 2.9 Gyr; Carretta et al. 2000), 
for Type la supernovae (14.9 ± 1.4 Gyr; Perlmut- 
ter et al. 1999), and from the shape of acoustic 
peaks in the cosmic microwave background spec- 
trum (14.0 ± 0.5 Gyr; Knox, Christcnscn, & Sko- 
rdis 2001). 

One advantage of the thorium chronology is 
that, once the initial and current values of Th/r 
in the star are specified, the age of the star de- 
pends only on the half- life of 232 Th determined in 
the laboratory. That is, one is not forced to in- 
voke complex (and no doubt incomplete) models 
of Galactic chemical evolution, although insight 
from making simple attempts can still be obtained 
(e.g., Cowan et al. 1997). Here, r denotes one of 
the stable r-process elements, often taken to be 

1 For convenience of nomenclature, in the context of radioac- 
tive cosmochronology it is common practice to refer to the 
"age" of a given star as the time difference between the 
production event (or events) of the radioactive (and sta- 
ble r-process) species that have been incorporated into the 
outer atmosphere of a low-mass (< 0.8Mq) star that has 
not yet evolved past the giant branch stage of evolution, 
and hence is amenable to spectroscopic study. 



Eu. In the most simple model, the initial value of 
Th/r may be regarded as that arising from a single 
supernova that exploded in the near vicinity (and 
perhaps may have even triggered the formation of) 
the presently observed star, e.g., CS 22892-052, 
early in the history of the Galaxy. In fact, the 
large dispersion of Eu/Fe observed in halo stars 
(more than 2 orders of magnitude) has been natu- 
rally explained by chemical evolution models that 
make such assumptions (Ishimaru & Wanajo 1999; 
Tsujimoto, Shigeyama, & Yoshii 2000). 

Thus far, the initial production of Th/r has 
been determined by fitting theoretical nucleosyn- 
thesis results to the solar r-process pattern, with 
the assumption that the r-pattern was univer- 
sal in all astrophysical environments (Cowan et 
al. 1999; Goriely & Clerbaux 1999). However, 
a recent abundance analysis of CS 22892-052 
with Keck I (HIRES) has shown that the lighter 
neutron-capture elements (Z < 48) are distinctly 
underabundant with respect to the solar pattern, 
which is well fit by species with Z > 56 (Sneden 
et al. 2000). This implies that there likely exists 
at least two r-process sites that must be invoked 
to account for the entire set of elements associ- 
ated with the r-process (Ishimaru & Wanajo 2000; 
Qian 2000). Furthermore, it should be noted that 
although Johnson & Bolte (2001) have derived 
an estimate of the minimum age of the universe, 
11.0 ± 4.2 Gyr, based on the abundances of tho- 
rium and other stable neutron-capture elements 
in five metal-poor stars, the estimated ages of 
these stars are very widely spread, from w 6 to 
18 Gyr (with overlapping error bars). In addition, 
recent spectroscopic studies of several extremely 
metal-poor halo stars with SUBARU (HDS) has 
revealed that the spread of estimated chronomet- 
ric ages from the Th-Eu pair appears to exceed 
the abundance uncertainties (Honda et al. 2002, 
in preparation). Taken at face value, these results 
imply that the universality of the r-pattern may 
not hold over the actinides either. 

The second discovered r-process-enhanced, ex- 
tremely metal-poor star, CS 31082-001 ([Fe/H] = 
—2.9) has provided a new, potentially quite pow- 
erful cosmochronometer, uranium (Cayrel et al. 
2001; Hill et al. 2002). In principle, uranium 
might be expected to be a more precise chronome- 
ter than thorium owing to its shorter half life. 
Cayrel et al. (2001) concluded that the age of 
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this star is 12.5 ± 3 Gyr, based on comparisons 
of the observed ratios of U/Th, U/Os, and U/Ir 
with their estimated production ratios (Cowan et 
al. 1999; Goricly & Clcrbaux 1999). The discov- 
ery of uranium in CS 31082-001 has also prompted 
additional laboratory work to refine the oscillator 
strengths, for both uranium and thorium (Nilsson 
et al. 2002a, b), that are required to turn spec- 
troscopic observations into measured abundances. 
However, analysis of this star has also provided 
conclusive evidence that the r-process pattern is 
not universal over the actinidcs; the ratio Th/Eu 
in this star is higher than that of the solar re- 
process ratio. CS 31082-001 would be younger 
than the solar system if the initial Th/Eu were 
taken to be universal, in conflict with the ages 
derived from ratios involving U and other sta- 
ble r-process elements close to the third r-process 
peak. Therefore, any age estimates that demand 
assumption of the universality of the r-process 
pattern may in fact be unreliable. 

In addition to the above non-universality prob- 
lem, the initial r-process pattern has thus far 
been determined theoretically by the superposi- 
tion of nucleosynthesis results, where one is forced 
to assume constant temperatures, neutron num- 
ber densities, and exposure times (or the num- 
ber of captured neutrons) - the most recent ex- 
ploration of this approach is presented in Schatz 
et al. (2002). These approximations have been 
necessary because of the lack of a reliable astro- 
physical model for the r-process site (at least that 
associated with the heavy r-process elements) . At 
the moment, the neutrino- wind scenario is consid- 
ered to be the most promising, ever since the work 
by Woosley et al. (1994), although alternative 
models involving neutron-star mergers have also 
shown some promise (Freiburghaus, Rosswog, & 
Thielemann 1999b). However, even the neutrino- 
wind model encounters serious problems, e.g., the 
high entropy (~ 400^^) that led to the robust r- 
processing in the Woosley et al. scenario has not 
been duplicated by other work (Takahashi et al. 
1994; Qian & Woosley 1996). 

The purpose of this paper is, therefore, to con- 
struct a realistic r-process model based on the 
neutrino-wind scenario, in order to derive initial 
production ratios that might be useful for estima- 
tion of the minimum age of the universe, in par- 
ticular by use of the U-Th chronometer pair. Re- 



cently, Wanajo et al. (2001) have demonstrated 
that the solar r-process pattern for A w 120 — 200 
is well reproduced in such a model from a compact 
proto-ncutron star, owing to the inclusion of gen- 
eral relativistic effects. Although recent hydrody- 
namic studies have encountered difficulties attain- 
ing the required compactness, we use their model 
along with some improvements to be discussed in 
§ 2. Although some severe problems yet remain, 
it is clearly important to determine the initial r- 
process pattern based on site-specific, rather than 
only site- independent, models. In § 3 the r-process 
nucleosynthesis calculation is described, and the 
results are compared to the solar r-process pat- 
tern. The initial electron fraction, Y e , is varied 
from 0.39 to 0.49, but taken constant with time 
during the neutron-capture process. In § 4, the 
age of CS 31082-001 is derived by comparison of 
the nucleosynthesis results obtained for each Y e 
with the abundance pattern of heavy elements in 
this star. The implications of this study are dis- 
cussed in §5. 

2. NEUTRINO WINDS 

Neutrino winds, in which the free nucleons ac- 
celerated by the intense neutrino flux near the neu- 
trino sphere of a core-collapse supernovae assem- 
ble to heavier nuclei, has been believed to be the 
most promising astrophysical site of the r-process 
since the work by Woosley et al. (1994). The 
rather high entropy (~ 400 A^fc) of the ejected 
matter in their hydrodynamic simulations led to 
production of the r-process elements that is in re- 
markable agreement with the solar r-process pat- 
tern. However, subsequent analysis of both hy- 
drodynamic and analytic approaches suggests that 
such high entropy may not be achieved in neutrino 
winds (~ 100 N A k; Takahashi et al. 1994; Qian & 
Woosley 1996). In addition, an unacceptable over- 
production of the elements with A ss 90 appeared 
in the r-process yields of Woosley et al. (1994). 

Recently, Otsuki et al. (2000) have studied the 
physical conditions required for the r-process in 
detail, using a semi-analytic model of a spheri- 
cal, steady neutrino wind, taking general relativis- 
tic effects into account (see also Qian & Woosley 
1996; Cardall & Fuller 1997; Thompson, Bur- 
rows, & Meyer 2001). They suggested that ro- 
bust physical conditions for the r-process were ob- 
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tained only if the proto-neutron star was as mas- 
sive as <~ 2.0Mq. Furthermore, Wanajo et al. 
(2001) showed from nucleosynthesis calculations 
that the r-process yields for the compact proto- 
neutron star models (with a mass of 1.9 — 2.0M Q 
and radius of 10 km) were in good agreement 
with the solar r-process pattern. This agreement 
was achieved due to the short dynamical timescale 
(< 10 msec) and moderately high entropy (> 140 
N^k) obtained at the early phase of the neutrino 
wind, which prevents most of the free nucleons to 
assemble into heavy nuclei. There are, however, 
some severe problems in their models. For exam- 
ple, it would appear difficult to obtain such a com- 
pact proto-neutron star with hydrodynamic mod- 
els of core-collapse supernovae using the equations 
of state of high-density matter that are currently 
available. This model also suffered from significant 
overproduction of elements with A w 90. Hence, 
although problems remain, we have chosen to be- 
gin with these models, and attempt to refine them 
in order to see if progress can be made toward the 
goal of obtaining more realistic production ratios 
of heavy elements that might prove useful for the 
U- and Th- chronologies. 

In this paper we use a model with a proto- 
neutron star with (gravitational) mass of 2.0M Q 
and radius of 10 km, as in Wanajo et al. (2001), 
whose nucleosynthesis results are in good agree- 
ment with the solar r-pattern for nuclei A « 
120 — 200. We briefly describe the models here, 
and point out some improvements that have been 
added. 

The system is treated as time stationary and 
spherically symmetric, and the radius of the neu- 
tron star is assumed to be the same as that of 
the neutrino sphere. The physical variables in 
the neutrino wind are then functions of the ra- 
dius only. The ejected mass by neutrino heat- 
ing is assumed to be negligible compared to the 
mass of the neutron star. Therefore, the gravi- 
tational field in which the neutrino-heated mat- 
ter moves can be treated as a fixed-background 
Schwarzschild metric. Time variations of the tem- 
perature and density are then solved by use of re- 
lations based on baryon, momentum, and energy 
conservation (equations (l)-(3) in Wanajo et al. 
2001). The source term in the equation of energy 
conservation is due to both heating and cooling by 
neutrino interactions. The gravitational redshift 



of the neutrino energies, and the bending of the 
neutrino trajectories expected from general rela- 
tivistic effects, are explicitly taken into account. 

For the cooling by annihilation of electron- 
positron pairs into neutrino-antineutrino pairs, 
the accurate table of Itoh et al. (1996) is uti- 
lized. All other neutrino heating and cooling rates 
are taken from Otsuki et al. (2000). The neu- 
trino luminosities, L u , of all neutrino flavors are 
assumed to be equal, and the RMS average neu- 
trino energies are taken to be 12, 22, and 34 MeV, 
for electron, anti-electron, and the other flavors of 
neutrinos, respectively. The equation of state for 
the electron and positron gas includes arbitrary 
relativistic pairs, which are of importance during 
en-processing. The electron fraction is set to be 
Y e « 0.43, which was mainly determined by the 
equilibrium of neutrino capture on free nucleons 
with the above luminosities and energies of the 
electron and anti-electron neutrinos at the neu- 
trino sphere. 

In this paper a freezeout temperature, Tf, is in- 
troduced to specify the final temperature of the 
neutrino winds. As can be seen in realistic hydro- 
dynamic simulations of neutrino winds, all trajec- 
tories specific to the r-process asymptotically ap- 
proach the same temperature in the stalled shock 
(Woosley et al. 1994). As discussed in § 3, this 
assumption is also reasonable from the nucleosyn- 
thctic point of view. Otsuki et al. (2000) set the 
mass ejection rate, M, at the neutrino sphere of 
each wind so that the temperature at 10000 km 
in radius is 0.1 MeV. This assumption may be in- 
adequate, as the physical conditions of the neu- 
trino sphere and the outer boundary are not nec- 
essarily causally connected. On the other hand, 
Wanajo et al. (2001) set the mass ejection rate to 
be M — 0.995 x M c , where M c is the maximum 
mass ejection rate at the neutrino sphere appropri- 
ate for the local physical conditions (the wind with 
M c becomes supersonic through the sonic point, 
while those with M < M c are subsonic through- 
out). However, the asymptotic temperatures dif- 
fered significantly from wind to wind, which is un- 
likely in a realistic situation. In this paper, there- 
fore, we have specified both M and Tj to avoid 
the above problems. The mass ejection rate is set 
to be M = M c . In fact, recent hydrodynamic sim- 
ulations show that the winds become supersonic, 
and then, in turn, decelerate to become subsonic 
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from the effects of the reverse shock coming from 
the outer iron mantle (Burrows, Hayes, & Fryx- 
ell 1995). Tf is set to be the same for all winds. 
The temperature, T, of each wind is thus replaced 
with Tf for all T < T f . The density is also re- 
placed to that obtained at T — Tf. The value of 
Tf is determined in § 3. 

3. THE r-PROCESS 

The r-process nucleosynthesis calculation, adopt- 
ing the model described above for the physical 
conditions, is obtained by application of an exten- 
sive nuclear-reaction network code. The network 
consists of ~ 3600 species, all the way from single 
neutrons and protons up to the fermium isotopes 
(Z = 100, see Figure 2). We include all relevant 
reactions, i.e., (n,j), (p,-f), (a, 7), (p,n), (a,p), 
(a, n), and their inverses. Reaction rates are taken 
from Thielemann (1995, private communication) 
for nuclei with Z < 46 and from Cowan, Thiele- 
mann, & Truran (1991) for those with Z > 47. 
The latter used the mass formula by Hilf et al. 
(1976). 

Our network does not include the nuclei near 
the neutron drip line for Z < 10, since no reli- 
able data are yet available. Terasawa et al. (2001) 
demonstrated that r-processing was hindered sig- 
nificantly by reaction flows near the neutron drip 
line for Z < 10, owing to their rather low asymp- 
totic temperature (— 6.2 x 10 8 K). As we sug- 
gest below, however, it may be unlikely that the 
asymptotic temperature is, in reality, less than 
~ 8 x 10® K. Thus, inclusion of these reactions 
would not significantly change our conclusions. 

The three-body reaction a(an,-f) 9 Be, which is 
of special importance as the bottleneck reaction 
to heavier nuclei, is taken from the recent exper- 
imental data of Utsunomiya et al. (2001). The 
weak interactions, such as /3-decay, /3-delayed neu- 
tron emission (up to three neutrons), and electron 
capture are also included, although the latter is 
found to be unimportant. Note that neutrino- 
induced reactions are not included. In fact, neu- 
trino capture on free nucleons might be of impor- 
tance, as they may serve to increase the electron 
fraction during seed element production (Meyer et 
al. 1998), while those on heavy nuclei are likely to 
be unimportant. However, it would be meaning- 
less to consider these effects without taking the 



accurate treatment of neutrino transport near the 
neutrino sphere into account. 

The a-decay chains and spontaneous fission 
processes are taken into account only after the 
freezeout of all other reactions. For the latter, 
all nuclei with A > 256 are assumed to decay by 
spontaneous fission only. The few known nuclei 
undergoing spontaneous fission for A < 256 are 
also included, along with their branching ratios. 
Neutron-induced and /3-delayed fissions, as well as 
the contribution of fission fragments to the lighter 
nuclei, are neglected. Obviously, these treatments 
of the fission reactions are oversimplified. Never- 
theless, this may be acceptable, at least to first 
order, as we demonstrate below. We leave the 
more accurate treatment of these matters to fu- 
ture work. In particular, /3-delayed fission may be 
of importance, as shown by Schatz et al. (2002). 

Each calculation is started at T 9 = 9 (where 
T 9 = T/10 9 K). At this point the nuclear sta- 
tistical equilibrium (NSE) consists mostly of free 
nucleons. The initial mass fractions of neutrons 
and protons are thus given by X n = 1 — Y e and 
X p = Y e , respectively. In Wanajo et al. (2001), 
the initial electron fraction Y e was fixed to be 0.40. 
In this paper Y e is varied from 0.39 to 0.49 to ex- 
amine its effect on the r-process. Note that, for 
simplicity, for each Y e the same trajectories are 
used, being derived with Y e w 0.43. The small 
change of Y e is, however, found to be unimportant 
for the temperature and density histories. 

The freezeout temperature, Tf, is determined 
as follows. At present, there is no positive evi- 
dence that the r-process pattern between the r- 
process-enhanced stars (e.g., CS 31082-001) and 
the solar abundances differs for Z > 56, except 
for Th and U (Hill et al. 2002). Thus, since our 
main goal is to estimate reliable yields of U and 
Th, the local r-process pattern near the third peak 
( 195 Pt) is treated most carefully, in order to repro- 
duce the solar r-process pattern. In Figure 1 the 
nucleosynthesis results before and after a-decay 
and fission (thick line), obtained with different 
Tgf(= Xy/10 9 K), are compared to the solar r- 
process pattern (points; Kappeler et al. 1989) nor- 
malized to the height of the third-peak nuclei. L v 
and Y e are taken to be 10 518 ergs s _1 and 0.43, 
respectively. The thin line denotes the results ob- 
tained by turning off the process of /3-delayed neu- 
tron emission. Note that the local pattern of the 
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r-element abundances (i.e., near A = 195) is not 
significantly altered as the result of changing L v 
and Y e (see Figures 4 and 5). Thus, the discus- 
sion below may hold even for superposition of the 
yields from multiple trajectories. 

For T 9 f = 0.6 - 0.8 (Figures la-c), the posi- 
tion of the third peak shifts to the low side by 
two units from A = 195. Moreover, the valley on 
the lower side of the third peak seen in the solar re- 
process pattern (A sa 183) is substantially filled up 
with the jagged pattern. On the other hand, for 
Tgf = 1.2 — 1.4 (Figures lg-i), the third peak sig- 
nificantly shifts to the high side from A = 195. In 
addition, the abundance pattern is unacceptably 
jagged compared to the solar r-process pattern. In 
this regard, the results with Tgf — 0.9 — 1.1 repro- 
duce the properties of the solar r-pattern near the 
third peak quite well (Figures ld-f). Indeed, the 
third peak is correctly located at A — 195. The 
pattern in the valley at A ss 183 is also reasonably 
reproduced. In this paper, therefore, we set Tgf 
to be 1.0, so that the predicted r-process pattern 
is in the best agreement with that of the solar re- 
process pattern near the third peak. Note that 
the second peak (A — 130) nuclei in Figure 1 are 
substantially underproduced compared to the so- 
lar r-abundances in this specific trajectory, which 
originates mainly from the higher L v trajectories 
(Figures 4 and 5). 

In Table 1, the abundance ratios among Eu, Ir, 
Pb, Bi, Th, and U are listed along with the isotope 
number fraction of 151 Eu. Note that, in this table, 
"Pb" means the sum of 206 Pb, 207 Pb, 208 Pb, and 
235 U, while "U" means 238 U, owing to the short 
mean half-life of 235 U compared to the cosmic age. 
As can be seen, the ratio U/Th is insensitive to the 
change of Tgf, except for the higher temperature 
cases with Tgf > 1.2. This is simply due to the 
fact that Th and U are close to one another in 
mass number. Thus, small variations of Tgf from 
1.0 may not change the chronometric ages sub- 
stantially, as long as the U-Th pair is employed. 
In contrast, other abundance ratios are found to 
be quite sensitive to the change of Tgf. For exam- 
ple, the ratio Th/Eu varies by more than a factor 
of two in the same temperature range. This im- 
plies that the chronometer pairs Th-r may be un- 
reliable unless the astrophysical model is strictly 
constrained. It is interesting to note that the ra- 
tio between Pb-Bi and Th-U can be changed by a 



factor of two even though they are mainly a-decay 
products. 

The isotope ratio of Eu can provide another 
constraint on the determination of Tgf. Recent 
spectroscopic studies have shown that the isotope 
number fractions of a few extremely metal-poor 
stars, including CS 22892-052 (Sneden et al. 2002) 
and CS 31082-001 (Aoki et al., in preparation) 
are close to the solar value of 0.478. This is also 
in good agreement with the calculated values for 
Tgf < 1.0 in the last column of Table 1. 

The strong influence of Tf on the final predicted 
r-process pattern can be recognized as the con- 
sequence of the frcczcout effects. Figure 2 shows 
the post-freezeout evolution of the r-process abun- 
dances for T 9 f — 1.0 with L v = 10 518 ergs s _1 and 
Y e = 0.43. The abundances are color-coded with 
their pattern as a function of mass number in the 
upper left. The dots denote the nuclei included 
in the reaction network, with the stable and long- 
lived nuclei shown as larger dots. At the freezeout 
of r-processing, which occurs 0.48 seconds after 
Tg = 9, a clear odd-even effect appears, owing to 
the (n,7) — (7,n) equilibrium (Figure 2a). After 
0.58 seconds, the abundances near the third peak 
are still along the neutron magic number N = 126 
(Figure 2b). As a result, the third peak slightly 
shifts to the larger side in mass number, due to 
the small fraction of neutrons emitted by photo- 
dissociation of fragile nuclei (e.g., A ~ 183) and 
by the /3-delayed process; these are quickly ab- 
sorbed by other nuclei (including those near the 
third peak). The abundance valley at A « 183 is 
then sharpened. At the same time the abundance 
pattern is significantly smoothed. At some point 
the abundances near the third peak separate from 
N = 126 due to /3-decay, becoming more and more 
robust against photo-dissociation, as well as with 
the decrease in the /3-delayed process (Figure 2c). 
The position of the third peak is thus fixed. The 
abundance pattern is smoothed further by the con- 
tinuation of non-equilibrium neutron emission and 
absorption among fragile nuclei. It is important 
to note that the separation from the N = 126 line 
takes place at earlier times for lower Tf, and vice 
versa, which determines the location of the third 
peak, as well as the smoothness of the pattern near 
it. Although the /3-delayed process is found to be 
unimportant for setting the locations of the sec- 
ond (A w 130) and the third (A w 195) peaks, it 
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does contribute to the enhancement of nuclei with 
A > 220, in particular for Tg/ > 0.9, as can be 
seen in Figure 1 . This results in the enhancement 
of Th-U/Pb-Bi for T 9/ - 1.0 found in Tabic 1. 

In order to calculate the mass-integrated re- 
process yields, the time evolution of the neu- 
trino luminosity, L v , is assumed to be L v (t) = 
Luoit/to)" 1 ■ L u0 and t are taken to be 10 52 6 ergs s 
(~ 4x 10 52 ergs s _1 ) and 0.2 s, respectively, being 
in reasonable agreement with the hydrodynamic 
results of Woosley et al. (1994). The time evolu- 
tion of the neutrino sphere radius is also assumed 
to be R v (t) = R v0 (t/t )~ 1/3 , where R v0 = 15 km. 
R v is replaced with 10 km for R v < 10 km. This 
rather rapid contraction of R v is required for ro- 
bust r-processing to occur at early times with 
a compact proto-neutron star model (Wanajo et 
al. 2001; Thompson, Burrows, & Meyer 2001). 
Figure 3 shows R v (km) and M c (M Q s^ 1 ), as 
well as the entropy per baryon S (AT^fc) and the 
timescale of seed abundance production r (msec) , 
as functions of L v . The timescale, r, is defined as 
the duration from 6 to 2.5 in Tg, which approxi- 
mately corresponds to the phase of a-processing, 
i.e., the seed abundance production prior to re- 
processing. Note that our definition of r differs 
from the e-folding time of the density often taken 
in other work, since the latter is not necessar- 
ily a good indicator of the seed abundance pro- 
duction. For L v « 5 — 10 ergs s _1 , r is signif- 
icantly short (« 5 — 10 msec) and S is moder- 
ately high (« 140 — 160A^4fc), which are suitable 
conditions for the r-process (preventing consump- 
tion of free neutrons prior to r-processing). For 
L v > 10 52 ergs s _1 , r increases and S decreases 
rapidly with L v , owing to the increasing R v . The 
total ejected mass, calculated by following the 
time evolution of L v and M c , is 5.0 x 10~ 4 M Q , 
more than one order of magnitude smaller than 
that of Woosley et al. (1994). This is perhaps due 
to the fact that our assumption of the steady-state 
flow docs not hold well for L v > 10 52 ergs s _1 . 
This high luminosity phase is, however, unimpor- 
tant for the r-process (Figure 4). 

The nucleosynthesis results for different L v 
with Ye = 0.43, before a-decay and spon- 
taneous fission, are shown in Figure 4. For 
logL„(ergs s -1 ) = 52.6, 52.4, and 52.2, no sig- 
nificant r-processing beyond the second peak oc- 
curs (Figures 4a-c). This is due to the low S 



(= 61, 81, HOiVyifc), regardless of the short r (= 
8.9, 7.3, 5.9 msec). For logT„(ergs s^ 1 ) = 52.0, 
51.8, and 51.6, robust r-processing occurs, pro- 
ducing the second- and third-peak nuclei, and 
substantial actinides (Figure 4d-f). This is a con- 
sequence of the rather short timescale of seed 
abundance production (r = 5.5, 7.5, 11 msec) 
and the high entropy per baryon (S = 141, 154, 
and 169 N A k). For iogT„(ergs s" 1 ) = 51.4, 51.2, 
and 51.0, the second and the third peaks are still 
prominent, but production of the actinides is not 
of importance with high S (= 184, 200, 217 N A k) 
and long r (= 17, 26, and 41 msec, Figure 4g- 
h). Note that, for logL^(ergs s" 1 ) = 52.0 - 51.0, 
the positions of the second (A w 130) and third 
(A « 195) peaks and the valley at A « 183 do not 
vary significantly (nor do the local patterns close 
to these features), although the global patterns 
are different (Figures 4d-i). This is due to the use 
of the same Tf (= 1.0 x 10 9 K) in all situations. 

We now integrate the nucleosynthesis results of 
trajectories with logL„(ergs s _1 ) = 52.6 — 51.0, 
using the time evolution of L v with M c (Figure 3). 
The interval of L v is taken to be 0.1 dex. These L v 
correspond to times t — 0.2 — 8.0 sec at the neu- 
trino sphere. The abundances of species formed 
after this period are neglected because of the small 
M c (< 2 x 10~ 6 M Q s" 1 ). In Figure 5, the mass- 
integrated abundances are compared with the so- 
lar r-process pattern (Kappeler et al. 1989) scaled 
to match 153 Eu. The a-decay and spontaneous 
fission after freezeout of all other reactions are in- 
cluded. The initial electron fraction, Y e , is varied 
from 0.39 to 0.49 to demonstrate its effect on the re- 
process. Note that in our models Y e is taken to be 
constant with time. Although Y e should certainly 
be a function of time, it is meaningless to con- 
sider this refinement without taking into account 
accurate treatment of neutrino transport near the 
neutrino sphere. As can be seen in Figure 5, for 
Y e = 0.39 — 0.43, the calculated abundances are 
in reasonable agreement with the solar r-process 
pattern between A w 120 and 200. For larger 
neutron excess, with Y e = 0.39 - 0.41, 232 Th and 
235,238tj are significantly overproduced compared 
to the solar r-process pattern (Figures 5a-c), as 
are other (stable) nuclei near the third peak. For 
Y e = 0.42, overall agreement of the abundance 
pattern oi A > 120 to the solar r-process pat- 
tern can be seen (Figure 5d). For Y e = 0.43, Th 
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and U are significantly underproduced compared 
to the solar r-process pattern (Figure 5e). In the 
current models, therefore, Y e should be < 0.42 to 
achieve sufficient production of Th and U. Note 
that the location of the third (A w 195) peak and 
the valley at A « 183 still hold even after inte- 
gration of the yields. The deficiency of nuclei at 
A = 130 — 140, seen for all Y e cases, could be due 
to our insufficient treatment of the physics of the 
neutrino winds for L v > 10 52 ergs s -1 , in which 
the assumption of the steady state winds is not 
manifest at all. More likely, this might be due 
to the nuclear properties of the mass model (Hilf 
ct al. 1976) employed in this study, as shown by 
Freiburghaus et al. (1999a). In their paper, calcu- 
lations with other mass formulae, FRDM (Mollcr 
ct al. 1995) and ETFSI (Aboussir ct al. 1995), were 
able to well reproduce the solar r-abundances in 
this mass region. 

As can be seen in Figure 5, overproduction of 
nuclei with A « 90 — 110 is evident, except for 
the Y e = 0.49 case. This is a consequence of the 
large M c for L v > 10 52 ergs s _1 (Figure 3), in 
which little r-processing occurs (Figures 4a-c) . For 
high Y e (= 0.49, Figure 5f), the overproduction 
significantly diminishes. This arises because, in 
the nuclide chart (see Figure 2), F eiSce( j > 0.446 
lines have intersections with a separation ener- 
gies of w 6 MeV on the lower side of the neu- 
tron magic number N = 50 (Hoffman et al. 1996; 
Freiburghaus et al. 1999a). At the same time, 
however, the r-process abundance of A « 120 — 200 
decreases substantially. This is due to the small 
levels of the seed abundance, owing to the ex- 
tremely strong a-rich freezeout by the less effi- 
cient a(an,j) 9 Be reaction with low neutron ex- 
cess (Hoffman et al. 1997). 

In Table 2, the total ejected masses of r-process 
abundances (A > 120) are shown, along with 
that of Eu. For Y e = 0.49, the ejected mass is 
about one order of magnitude smaller than for 
other Y e cases. According to the chemical evo- 
lution study of the Galactic halo by Ishimaru & 
Wanajo (1999), ~ 5 x 10~ 7 M© of Eu per r-process 
event is needed to reproduce the large dispersion 
of observed Eu/Fc in extremely metal-poor stars. 
The mass of Eu obtained for Y e — 0.49 is too small 
compared to that level, while remaining in reason- 
able agreement for lower Y e cases. We speculate 
that, if the primary r-process site is really neutrino 



winds, as discussed in this study, the electron frac- 
tion changes from Y e « 0.49 to < 0.42 rapidly with 
time. In this case, the overproduction would not 
appear at the early phase [L v > 10 52 ergs s _1 ), 
but would be followed by later r-nuclei (A > 120) 
production (L v < 10 52 ergs s _1 ). As can be seen in 
the third column in Table 2, except for Y e = 0.49, 
the isotopic fraction of 151 Eu is in good agreement 
with both the solar value and the observational es- 
timates (Sneden et al. 2002; Aoki et al., in prepa- 
ration), even after mass-integration of the yields. 

The last column of Table 2 shows the mass of 
fission fragments, whose contribution to the lighter 
nuclei is neglected in this study. This assump- 
tion appears adequate since the contribution of 
fission fragments is less than 5% of the total mass 
of nuclei with A > 120. This follows since, in 
this study, the neutron-to-seed ratio at the onset 
of r-processing (Tg = 2.5) is 138 at most, with 
(A^eed) ~ 95. This results in no significant pro- 
duction of nuclei with A > 255, where fission dom- 
inates. 

4. THE U-Th COSMOCHRONOLOGY 

In Figure 6 the available spectroscopic abun- 
dance data for CS 31082-001 (dots) are compared 
with the nucleosynthesis results discussed in § 3 
(thick line) and with the solar r-pattern (thin 
line), scaled at Eu (Z = 63). The data for the 
neutron-capture elements in this star are taken 
from Hill et al. (2002). CS 31082-001 is the 
first extremely metal-poor star ([Fe/H] = —2.9) in 
which the UII line (3859.59 A) has been detected 
(Cayrel et al. 2001). The abundance pattern of 
this star between Ba (Z = 56) and Ir (Z = 77) 
is in excellent agreement with the solar r-process 
pattern, except for a few elements (Tb, Hf, Os). 
It is possible that, in these few cases, the disagree- 
ments are due to difficulties in the abundance es- 
timates (Hill et al. 2002). Note that some lighter 
elements (Rh, Ag) are not in agreement with the 
solar r-process pattern, as has also been seen in 
CS 22892-052 (Sneden ct al. 2000), and which im- 
plies the presence of two distinct r-process sites 
(Ishimaru & Wanajo 2000; Qian 2000). 

Nucleosynthesis results for all Y e cases are in 
good agreement with the abundance pattern of 
CS 31082-001 for the elements between Pr (Z = 
59) and Tm (Z = 69), with the exception of Tb. 
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Note that a few of the predicted abundances for el- 
ements near the second peak (Ba, La, and Ce) are 
significantly deficient compared to both the obser- 
vational and solar patterns. As mentioned in § 3, 
this might be due to the properties of the nuclear 
mass model employed in this study. In future work 
we plan to check if these elements are reproduced 
with other mass formulae, since the observational 
data suggest the universality of the r-process pat- 
tern between the second and the third r-process 
peaks. The platinum-peak elements Os and Ir are 
in reasonable agreement, except for the Y e = 0.49 
case, although the height of the 195 Pt differs. In 
addition, the large abundance of Os in CS 31082- 
001 compared to the scaled solar r-process pattern 
makes it questionable if the third peak is really lo- 
cated at Pt. It would clearly be useful to obtain 
measurements of the Pt abundance in this star, 
so that the height, as well as the position, of the 
third peak could be better constrained. Note that 
235 U is assumed to have a-decayed away because 
of its relatively short half life (0.704 Gyr). 

As discussed by Hill et al. (2002), the age of 
CS 31082-001 can be inferred by application of 
one or more of the following three relations: 

t* TKr = 46.67 pog(Th/r) - log(Th/r) now ] Gyl) 
t^ r = 14.84 [log(U/r) - log(U/r) now ] Gyr (2) 
*u,Th = 21.76 [log(U/Th) - log(U/Th) now ] 

with the half lives of 232 Th (14.05 Gyr) and 238 U 
(4.468 Gyr), r is a stable r-element, and the sub- 
scripts "0" and "now" denote the initial and cur- 
rent values derived by theory and observation, re- 
spectively In principle, r and Th can be 
more precise chronometers than t^ h r , owing to the 
smaller coefficients in front of equations (2) and 
(3). In practice, the U-Th pair (equation (3)) is ex- 
pected to be the most precise chronometer owing 
to the similar atomic properties of these species, 
which serve to diminish the observational errors 
associated with determination of their ratio. Fur- 
thermore, the ratio U/Th is less dependant on the 
theoretical model used, since these species are sep- 
arated by only two units in atomic number, as well 
as being mostly the result of a-decayed products 
(Table 1). As stable elements r, Eu, Os, and Ir are 
selected. The Eu abundance derived by spectro- 
scopic analyses is more reliable than for any other 



elements with Z > 59. The Pt-peak elements Os 
and Ir are co-produced with Th and U in the same 
(or close in L v ) trajectories, and thus are perhaps 
less model dependent as well. Note that, if possi- 
ble, the Pb abundance in this star should be used, 
since it is mainly synthesized by the same a-decay 
chains as Th and U. However, Pb has not yet been 
detected in CS 31082-001. and for now only an up- 
per limit to its abundance is available. For reasons 
that are not yet clear, the upper bound on Pb re- 
ported by Hill et al. (2002) is lower than the solar 
value, being satisfied only by the Y e < 0.42 cases 
in which Th and U are not significantly produced. 
Thus, Pb is omitted for the age determination in 
this study, although its future measurement would 
clearly be of importance. 

Table 3 lists the predicted initial abundance ra- 
tios among Eu, Os, Ir, Pt, Th, and U (A = 238) 
for all the Y e cases considered in the present study. 
The ratio U/Th is found to be insensitive to Y e 
(except for the extreme case of Y e = 0.49), which 
is in excellent agreement with those constrained 
by the solar r-process abundances in Goriely & 
Arnould (2001). It should be noted that these 
values can be applied not only to CS 31082-001, 
but also to other metal-poor halo stars in which 
uranium and thorium (and platinum) are detected 
by future observations (Cowan et al. 2002). With 
the use of these ratios as the initial values, the 
age of CS 31082-001 is derived as shown in Ta- 
ble 4, as well as in Figure 7. It is noteworthy 
that the age obtained from the U-Th chronometer 
pair is considerably more robust than the alter- 
natives, resulting in a range ~ 13.5 — 14.2 Gyr, 
with an exception of 12.2 Gyr for the Y e = 0.49 
case. Therefore, we use only the U-Th chronome- 
ter for dating of this star, while the others are 
regarded as consistency checks only. In particu- 
lar, the Th-Eu pair should be used to constrain 
the model parameter Y e in this study, the species 
being widely separated in mass number. For this 
reason, Y e is constrained to be « 0.40, owing to 
the consistency between ages based on the U-Th 
pair and others (especially Th-Eu). The ages ob- 
tained by use of Os are significantly higher than 
those by others, because of the previously noted 
discrepancies between observation and the theo- 
retical prediction (Figure 6). As seen in Table 4, 
the difference between the ages obtained from U- 
Th between Y e = 0.39 and 0.41 is only 0.1 Gyr. 
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On the other hand, the observational error asso- 
ciated with the U/Th ratio, 0.11 dex, leads to an 
uncertainty of 2.4 Gyr from application of equa- 
tion (3). We conclude, therefore, that the age of 
CS 31082-001 is 14.1 ± 2.5 Gyr. 

It should be emphasized that the age determi- 
nation made in the present paper is derived from 
the exploration of a single specific model based 
on the neutrino-wind scenario, and relies as well 
on the choice of one of the many existing nuclear 
data sets. Thus, it is clear that our final result 
may be changed if we were to make use of other 
astrophysical models. The nuclear mass models 
can also significantly affect the age determination 
(Goriely & Arnould 2001). Obviously, more study 
with other astrophysical models, as well as with a 
number of different nuclear data sets, are needed 
to derive the final results. 

5. SUMMARY AND CONCLUSIONS 

We have constructed a robust model of the r- 
process based on the neutrino-wind scenario for 
application to the U-Th cosmochronology. The 
model of neutrino winds was taken from Otsuki 
ct al. (2000) and Wanajo et al. (2001), along 
with some improvements described above. The r- 
process with various neutron star masses (Mns) 
and radii (w R v ) has already been explored ex- 
tensively by Wanajo et al. (2001). Thus, M NS 
and (asymptotic) R v were fixed to be 2.0M Q and 
10 km, respectively, for which the most robust re- 
processing was realized. On the other hand, the 
initial electron fraction Y e (assumed constant with 
time) was taken to be a free parameter, rather 
than fixed to be 0.40 as in Wanajo et al. (2001). 

The model was tuned to fit the solar r-process 
pattern, since recent spectroscopic studies of ex- 
tremely metal-poor stars have indicated a univer- 
sal (solar r-like) abundance pattern, at least be- 
tween the second and third r-process peaks. At- 
tention was paid to reproduce the pattern near 
A « 195, being close to Th and U in mass num- 
ber. We introduced the freezeout temperature, 
Tf, as the most crucial parameter for construct- 
ing the model. Comparing the nucleosynthesis 
results with the solar r-process pattern, Tf was 
determined to be 1.0 x 10 9 K. The nucleosynthe- 
sis yields were mass-integrated assuming simple 
evolutions of L v and R v . The results between 



A « 120 and 200 for Y e = 0.39 - 0.43 models 
were in good agreement with the solar r-pattern, 
although the production of Th and U differed from 
model to model. 

The results of our nucleosynthesis predictions 
were compared to the spectroscopic abundance 
pattern of the r-process-enhanced metal-poor star 
CS 31082-001. Explorations of various combina- 
tions between the actinides and stable isotopes in- 
dicated that the Y e = 0.40 case was found to be 
the best models for the age determination in this 
study. As a result, the age of CS 31082-001 was 
determined to be 14.1 ± 2.5 Gyr by use of the U- 
Th chronometer, which can be regarded as a hard 
lower limit on the age of the universe. In fact, 
this age is in excellent agreement with the age of 
the universe derived by other dating techniques 
(e.g., globular clusters, Type la supernovae, and 
the cosmic-microwave background) , as well as with 
other stellar radioactive age estimates using the 
Th-Eu (Sneden et al. 1996; Cowan et al. 1997) and 
U-Th (Hill ct al. 2002; Schatz et al. 2002; Cowan 
et al. 2002) pairs. Such a strong constraint on 
the age of this star was only made possible by the 
use of U and Th in conjunction with one another. 
Without uranium, the chronometers involving Th 
alone presented only rather loose constraints. 

In this study the Pb abundance in CS 31082- 
001 was not available for age determination, since 
only an upper limit on its abundance has been 
presented to date. It should be noted that the 
predicted Pb abundance for Y e = 0.40 case is 
higher than this upper limit, which seems diffi- 
cult to solve (Schatz et al. 2002). It is noteworthy, 
however, that freezeout effects can increase the ra- 
tio Th-U/Pb-Bi (see § 3), which should be exam- 
ined further, as well as with other nuclear data 
sets. It remains possible that the solar r-process 
Pb abundance, inferred from the residual of the 
solar s-process components, is not correct (Arlan- 
dini et al. 1999). Clearly, re-investigation of the 
Pb abundance in CS 31082-001, as well expansion 
of the searches for additional r-process-enhanced 
metal-poor stars, in which both U and Th are mea- 
surable, are highly desirable in the near future. It 
is also of importance to obtain measurements of 
the Pt abundance in CS 31082-001 (which is only 
possible with the strong Pt absorption features ob- 
servable in the near UV from space) , so as to con- 
strain the height, as well as the position, of the 
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third peak (Cowan et al. 2002). 

We emphasize that the conclusions in this study 
should be considered as merely the first step to- 
ward a better understanding of the dating tech- 
nique involving the U-Th chronometer. In partic- 
ular, we note that the inferred age of CS 31082-001 
could be significantly changed if we had utilized a 
different nuclear mass model, as demonstrated by 
Griely and Arnould (2001). We expect, however, 
that a detailed comparison of the nucleosynthesis 
results with the solar r-process pattern, as in § 3, 
can distinguish better mass models. We suggest 
that the freezeout temperature, Tj (or the asymp- 
totic temperature in more realistic models), is a 
useful parameter, as the freezeout effects are of 
importance to the final r-process pattern. While 
the waiting-point approximation, often referred to 
as the "classical r-process model," has been an ef- 
fective (and thus far, the only) method used to 
test the nuclear mass models, we should keep in 
mind that non-equilibrium processes can lead to 
non-negligible effects, as discussed in § 3 (e.g., the 
positions of the r-process peaks, the smoothness of 
abundance pattern, and the ratio Th-U/Pb-Bi). 

It should also be noted that we would not in- 
sist that the neutrino-wind model explored in this 
study is necessarily the best (or only) astrophys- 
ical site for the r-process. In fact, the physical 
conditions in this model, namely, the rather mas- 
sive proto-neutron star (2.0M Q ) with rapid con- 
traction of R v to be 10 km, seem difficult to be 
achieved according to the current hydrodynamic 
studies. Clearly, more study is needed to seek the 
true r-process site. Nevertheless, it is encouraging 
that the age of the CS 31082-001 determined in 
this study is in excellent agreement with the age 
of the universe derived by other dating techniques. 
We hope our results provide some fresh insights 
to the future modeling of the true astrophysical 
r-process site. 
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Table 1 
Abundance Ratios 



Tgf 


Th/Eu 


Th/Ir 


Th/Pb a 


Th/Bi 


U/Eu a 


U/Ir a 


U/Pb a 


U/Bi a 


U/Th a 


/( 151 Eu) b 


0.60 


0.24 


0.02 


0.08 


0.21 


0.11 


0.01 


0.04 


0.10 


0.44 


0.47 


0.70 


0.36 


0.02 


0.09 


0.27 


0.17 


0.01 


0.04 


0.13 


0.47 


0.46 


0.80 


0.55 


0.02 


0.13 


0.36 


0.25 


0.01 


0.06 


0.17 


0.46 


0.45 


0.90 


0.46 


0.02 


0.14 


0.50 


0.23 


0.01 


0.07 


0.25 


0.50 


0.44 


1.00 


0.38 


0.02 


0.15 


0.45 


0.19 


0.01 


0.07 


0.22 


0.49 


0.52 


1.10 


0.41 


0.03 


0.14 


0.47 


0.18 


0.02 


0.06 


0.21 


0.44 


0.41 


1.20 


0.17 


0.06 


0.09 


0.46 


0.06 


0.02 


0.03 


0.16 


0.36 


0.40 


1.30 


0.05 


0.09 


0.05 


0.23 


0.01 


0.03 


0.02 


0.07 


0.32 


0.24 


1.40 


0.03 


0.15 


0.06 


0.27 


0.01 


0.03 


0.01 


0.05 


0.18 


0.19 



ap b = 206p b + 207p b + 208p b + 235tJ ; TJ = 23& V ^ ^ 

b isotope number fraction of 151 Eu 



Table 2 
Ejected Mass (M q ) 





A > 120 




Eu 




/(151 Eu )a 


fission b 


0.39 


1.5 x 10" 


4 


6.2 x 10" 


-7 


0.50 


0.042 


0.40 


1.2 x 10" 


4 


5.4 x 10" 


-7 


0.50 


0.025 


0.41 


9.5 x 10" 


5 


4.3 x 10" 


-7 


0.52 


0.014 


0.42 


7.4 x 10- 


5 


4.5 x 10" 


-7 


0.51 


0.005 


0.43 


5.6 x 10" 


5 


4.6 x 10" 


-7 


0.51 


0.002 


0.49 


4.5 x 10" 


6 


6.8 x 10" 


-8 


0.40 


0.000 



a isotope number fraction of 151 Eu 

b mass fraction of fission fragments (^4 > 120) 
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Table 3 
Abundance Ratio (Log Scale) 





Th/Eu 


Th/Os 


Th/Ir 


Th/Pt 


U/Eu 


U/Os 


U/Ir 


U/Pt 


U/Th 


0.39 


0.18 


-0.18 


-0.60 


-1.08 


-0.11 


-0.47 


-0.89 


-1.37 


-0.29 


0.40 


0.05 


-0.41 


-0.84 


-1.27 


-0.24 


-0.70 


-1.13 


-1.57 


-0.29 


0.41 


-0.11 


-0.68 


-1.10 


-1.48 


-0.40 


-0.97 


-1.39 


-1.76 


-0.29 


0.42 


-0.58 


-1.16 


-1.56 


-1.85 


-0.89 


-1.47 


-1.87 


-2.16 


-0.31 


0.43 


-0.92 


-1.43 


-1.82 


-2.06 


-1.24 


-1.75 


-2.13 


-2.38 


-0.32 


0.49 


-2.75 


-2.50 


-2.83 


-3.17 


-3.13 


-2.89 


-3.21 


-3.55 


-0.38 



Table 4 
Ages of CS 31082-001 (Gyr) 



Y e 


Th/Eu 


Th/Os 


Th/Ir 


U/Eu 


U/Os 


U/Ir 


U/Th 


0.39 


18.77 


57.52 


27.18 


15.63 


27.95 


18.30 


14.16 


0.40 


12.61 


46.73 


16.09 


13.62 


24.47 


14.73 


14.10 


0.41 


5.17 


34.01 


3.64 


11.32 


20.49 


10.83 


14.19 


0.42 


-16.90 


11.67 


-17.55 


3.97 


13.05 


3.76 


13.70 


0.43 


-32.54 


-0.84 


-29.64 


-1.12 


8.96 


-0.20 


13.53 


0.49 


-118.21 


-51.05 


-76.97 


-29.30 


-7.94 


-16.18 


12.16 
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Fig. 1. — Comparison of the nucleosynthesis yields (line) for Tf = (a) 0.6, (b) 0.7, (c) 0.8, (d) 0.9, (e) 1.0, 
(f) 1.1, (g) 1.2, (h) 1.3, and (i) 1.4, with the solar r-pattern (dots), scaled at the height of the third peak 
for logZ/„ (ergs s -1 ) = 51.8 with Y e = 0.43. The thick and thin lines represent situations with and without 
/3-delayed neutron emission, respectively. 
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N •» 

Fig. 2. — Post-freezout evolution of the r-process yields at t — (a) 0.48, (b) 0.58, and (c) 2.13 seconds for 
\ogL v (ergs s _1 ) = 51.8 and Y e = 0.43. The abundances are color-coded in the nuclide chart. The pattern as 
a function of mass number is shown in the upper left of the each panel. The nuclei included in the reaction 
network are denoted by dots, with the stable and lorrg^lived isotopes represented by large dots. 
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Fig. 3. — The entropy per baryon (S), timescale of the seed production (t), neutrino sphere radius (R v ), 
and mass ejection rate (M c ) as functions of the neutrino luminosity (L v ), in units denoted along the axes. 
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Fig. 4.— Calculated yields for \ogL v (ergs s" 1 ) = (a) 52.6, (b) 52.4, (c) 52.2, (d) 52.0, (e) 51.8, (f) 51.6, (g) 
51.4, (h) 51.2, and (i) 51.0, with Y e = 0.43 
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Fig. 5.— Comparison of the mass-integrated yields (line) for Y e = (a) 0.39, (b) 0.40, (c) 0.41, (d) 0.42, (e) 
0.43, and (f) 0.49 with the solar r-pattern (dots) scaled at 153 Eu, as functions of mass number. 
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Fig. 6. — Comparison of the mass-integrated yields (the thick line) for Y e — (a) 0.39, (b) 0.40, (c) 0.41, 
(d) 0.42, (e) 0.43, and (f) 0.49, scaled at Eu (Z = 63), with the abundance pattern of CS 31082-001 (filled 
circles, with observational error bars), as functions of atomic number. For Pb, the observed upper limit is 
shown by the open circle. The scaled solar r-pattern is shown by the thin line. 
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Fig. 7. — Ages of CS 31082-001 derived from various chronometer pairs. The robustness of the U-Th pair is 
clearly shown. The superiority of the U-r pairs compared to those of Th-r can also be seen. 
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